Spatiotemporal vortex solitons in hexagonal arrays of waveguides 
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By means of a systematic numerical analysis, we demonstrate that hexagonal lattices of parallel 
linearly-coupled waveguides, with the intrinsic cubic self-focusing nonlinearity, give rise to three 
species of stable semi-discrete complexes (which are continuous in the longitudinal direction), with 
embedded vorticity S: triangular modes with S = 1, hexagonal ones with S = 2, both centered 
around an empty central core, and compact triangles with S = 1, which do not not include the empty 
site. Collisions between stable triangular vortices are studied too. These waveguiding lattices can 
be realized in optics and BEC. 
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I. INTRODUCTION 

Lattice solitons are a topic of great interest to on- 
going studies of nonlinear dynamics in photonic media 
and BECs (Bose-Einstein condensates) [1 j. These local- 
ized modes are produced by the interplay of the intrin- 
sic nonlinearity of the medium with an effective periodic 
potential induced in it by permanent or virtual lattice 
patterns. In fact, the lattice may itself be a nonlinear 
structure if it is induced by a spatially periodic mod- 
ulation of the local nonlinearity [2]. In the limit of a 
deep periodic potential, the fundamental models of lat- 
tice media reduce to various versions of the discrete non- 
linear Schrodinger (DNLS) equation [3]. The realization 
of the one-dimensional (ID) DNLS model in arrayed op- 
tical waveguides was originally proposed in Ref . [4 . The 
same model was later applied to BECs loaded into deep 
optical-lattice potentials [5] (see Ref. [6J for a brief re- 
view). A physical realization of the DNLS model is also 
possible in the form of lattices of microcavities which 
serve as traps for polaritons [7]. 

Lattice solitons take the form of discrete solitons in 
terms of the DNLS equations, which correspond to quasi- 
discrete solitons in the respective experimental settings. 
Such solitons were created in a set of semiconductor 
waveguides built on top of a slab substrate [8 j, and also 
in arrays of optical fibers [9] . In addition to using perma- 
nent photonic structures, quasi-discrete solitons were also 
made in virtual waveguiding arrays, using the versatile 
technique of inducing interference lattices in photorefrac- 
tive crystals [TO] , The latter method was used to create 
the first examples of 2D quasi-discrete fundamental soli- 
tons [11 j, which was followed by the making of vortex 
solitons [12] , i.e., localized lattice excitations with embed- 
ded vorticity, that were predicted in Ref. [13 . Another 
significant contribution to this area was the creation of 
2D solitons in a bundle of fiber-like waveguides written 
in bulk silica [14]. Such arrays and bundles are created 



by means of tightly focused femtosecond laser pulses [15] . 

Following the analysis of the fundamental localized dis- 
crete vortices with topological charge S — 1 [13 , their 
higher-order counterparts, with S > 1, and multipole 
discrete solitons, such as quadrupoles, were predicted in 
Refs. [16] [17]. Many other objects were studied in this 
area, including supervortices (circular chains of compact 
vortices with an imprinted overall topological charge, 
which is independent of the vorticity of the individual ed- 
dies [17]), necklace-shaped patterns [18], discrete solitons 
in hexagonal and honeycomb lattices [19 , 20 , composite 
semidiscrete spatial solitons in arrays of waveguides with 
quadratic and cubic nonlinearities [21] , quasi-discrete 
topological solitons in photonic-crystal fibers [22 , etc. 
Nonstationary soliton effects were studied too. These in- 
clude the mobility of discrete solitons |23] [24] , collisions 
between traveling ones [24 , 25 , and the onset of the spa- 
tiotemporal collapse in self-focusing arrayed waveguides 
[26]. 

Most works on lattice solitons dealt with the spatial- 
domain settings. In particular, optically-induced lattices 
in photorefractive crystals do not makes it possible to 
observe the evolution in the temporal domain because of 
a very large response time in these materials. However, 
the spatiotemporal dynamics can be realized in waveg- 
uiding arrays written in bulk silica [15], where the spa- 
tially localized quasi-discrete patterns in the transverse 
plane can be combined with the temporal self-trapping in 
the longitudinal direction. Recently, the creation of the 
corresponding quasi-discrete "light bullets" was reported 
in this system [27] (for a review of spatiotemporal solitons 
in nonlinear optics and BEC, see Ref. [28). Previously, 
a number of manifestations of the spatiotemporal self- 
trapping in similar settings were studied theoretically, 
including the related modulational instability [29 , for- 
mation of "bullets" in fiber arrays [30] and photonic wires 
[3T] , and self-compression [32] and steering [33] of pulsed 
beams. Continuing the work in this direction, semi- 
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discrete spatiotemporal surface solitons were introduced, 
as surface modes, in models of semi- infinite waveguide 
arrays [34], and in a system with an interface between 
different arrays [35]. Also analyzed were spatiotemporal 
solitons in waveguide arrays with the quadratic nonlin- 
earity [36] . 

Once stable fundamental spatiotemporal soliton com- 
plexes in bundled arrays of waveguides are available, it 
is natural to seek for vortex solitons in the same set- 
ting. A systematic analysis of spatiotemporal vortices 
and quadrupoles in the model based on the square lat- 
tice of discrete waveguides was reported in Ref. [37] . 
A vast stability area was found for the solitary vortices 
with S = 1 and quadrupoles, which are built as rhom- 
buses, alias on-site-centered modes, with respect to the 
underlying lattice (the rhombus is built as a set of four 
"bright" cores, with a nearly "dark" one at the cen- 
ter). The stability region is much smaller for the off-site- 
centered modes of the "square" type, without an empty 
pivotal site in the middle (the reduced stability domain 
of square-shaped vortices and quadrupoles, in compari- 
son with their rhombic counterparts, in a generic feature 
of topological solitons in lattice media [38). All the spa- 
tiotemporal vortex solitons with S =2 were found to be 
unstable unstable in the same model. Further, collisions 
between stable vortices and quadrupoles (with identical 
or opposite topological charges), propagating along the 
bundle in opposite directions, were analyzed in Ref. [39] . 
Four different outcomes of the collisions were identified: 
rebound of slowly moving solitary vortices, fusion, split- 
ting, and quasi-elastic interactions between fast ones. 

Hexagonal lattices may be created by means of the 
same techniques which were used for the building the 
square-shaped structures. On the other hand, the change 
of the underlying geometry may essentially alter funda- 
mental properties of topological lattice solitons [3] [19]. In 
particular, it was predicted theoretically and conformed 
in an experiment that spatial solitons in the form of dou- 
ble vortices (with S = 2) in hexagonal lattices may be 
stable, while their unitary counterparts (with 5 = 1) are 
unstable. 

The objective of the present work is to study spa- 
tiotemporal vortex solitons in hexagonal lattices of dis- 
crete waveguides. The model is formulated in Section 2, 
and at the end of it we also briefly consider fundamen- 
tal solitons, driven by a temporally self-trapped pulse 
in a single waveguiding core. In Section 3, we demon- 
strate, also in a brief form, that a straightforward input 
in the form of a hexagon-shaped spatiotemporal vortex 
with S = 1 always leads to a decay. Nevertheless, three 
different species of stable spatiotemporal complexes with 
the embedded vorticity are revealed by a systematic nu- 
merical analysis. In Section 4, we demonstrate that a 
spatiotemporal input of a triangular shape generates self- 
trapped vortices in the form of triangles with an empty 
core in the middle. Further, in Section 5 it is shown 
that a hexagonally shaped input with S = 2 produces 
stable spatiotemporal hexagons with the same (double) 



topological charge. Finally, a modified (shifted) input 
ansatz gives rise to stable densely packed triangular vor- 
tices with 5 = 1, without an empty central core, as shown 
in Section 6. In addition to the study of these species 
of spatiotemporal vortex solitons, in Section 7 collisions 
between counterpropagating triangular ones are studied. 
The paper is concluded by Section 8. 

II. THE MODEL AND FUNDAMENTAL 
SOLITONS 

We consider the hexagonal array of nonlinear waveg- 
uides, with cells in the transverse lattice numbered as 
shown in Fig. [I] The transmission of waves in the ar- 
ray is described by the following system of coupled NLS 
equations, written in the scaled form, similar to that used 
in many earlier works [29] -[37], [27] : 

'ld z U rn n -\- [Urn — i^ji— x -|- Urn—l,n -\- %l rn n — \ -\- U rn n J r \ 

m,nj 

+ (l/2)<9 £ 2 W m , n + Um^ |^ m ,n| 2 = 0. (1) 

In terms of the optical setting, z and t are, respectively, 
the propagation distance and reduced time, assuming 
that each guiding core features the anomalous chromatic 
dispersion and cubic self-focusing, while — /i is the propa- 
gation constant of the localized solution to be sought for. 
In terms of the corresponding BEC model, Eqs. is 
a system of coupled discrete Gross-Pitaevskii equations 
[40] , with z and t playing the roles of the scaled time and 
axial coordinate, respectively, while \i is the chemical po- 
tential. 

Simulations of Eq. ([I]) were carried out in the Fourier 
domain, with the help of the standard fourth-order 
Runge-Kutta scheme, the nonlinear term being evalu- 
ated by means of the combination of inverse and direct 
fast Fourier transforms at each sub-step of the scheme. 
We used an 11 x 12 matrix in the plane of (m,n), 512 
points for variable t in the computation window of width 
At — 20, and the stepsize in the propagation direction 
dz = 5 x 10 -4 . The use of the Fourier transform implies 
periodic boundary conditions in t, which make sense if a 
characteristic temporal size of the localized objects will 
be essentially smaller than At = 20. As concerns the 
boundary conditions for the discrete coordinates m and 
n, the values of u m ^ n corresponding to the coordinates 
which fall outside of the computation box are replaced 
by zeros. 

Before proceeding to the search for complex spatiotem- 
poral vortical patterns, it makes sense to test the prop- 
agation of fundamental solitons, which are carried, es- 
sentially, by a temporal pulse in a single core. For this 
purpose, the simulations were initiated with obvious ini- 
tial conditions, 

^o,o [z = 0) = r] sech (rjt) , (2) 
r) = V2 (6+70, (3) 
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FIG. 1. (Color online) The setup and notation: We con- 
sider the hexagonal array of cylindrical waveguides, as shown 
in the upper left corner. The transverse distribution of the 
light intensity of the propagating waves in the guiding cores 
is displayed symbolically in the bottom left corner. Each core 
is assumed to be a single-mode waveguide, represented by 
wave function Um,n, with discrete coordinates (m,n) defined 
as shown in the figure. The map of integers m, n into the 
Cartesian coordinates in the transverse plane, (x, y), is per- 
formed as per Eq. (J5|, a being the width of the hexagonal 
cell. 



setting u m , n (z = 0) = at \m\ + \n\ ^ 0. The simula- 
tions were run in interval 1 < \i < 16 of values of the 
propagation constant. 

It has been concluded that input ([2| decays, under 
the action of the lattice diffraction, at \i < 6.6, and a 
stable fundamental soliton, concentrated in the central 
core, is formed at /i > 6.6. The temporal pulse which lies 
at the core of the so created fundamental soliton is not 
quite stationary, but rather features regular pulsations, 
as shown in Fig. [2] 

One may surmise that the oscillations of the funda- 
mental soliton could be a result of its interaction with 
the radiation background, which was generated by the 
input field in the course of self-trapping into the funda- 
mental soliton. To check this possibility, the background 
around the soliton was explicitly removed, at a particu- 
lar step of the simulations. Nevertheless, the oscillations 
remain virtually unaffected by the "cleaning", i.e., they 
seem to be a genuine feature of the dynamics of the soli- 
ton, possibly representing its intrinsic mode. 



III. THE HEXAGONAL INPUT: A 
TRANSITION TO INSTABILITY 



First, we attempted to create hexagonal vortical modes 
with 5=1, which seems a natural approach to the 
system based on a hexagonal lattice. To this end, we 
used the following input, based on an ansatz factorized 
in the longitudinal (temporal) and transverse (spatial) 
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FIG. 2. (Color online) (a) Oscillations of the fundamental 
soliton generated by input J2|, |3| with \i — 8. (b) The up- 
per and lower curves illustrate the oscillations by showing, 
respectively, the largest and smallest values of the soliton's 
amplitude, i.e., max 2 (maxt (|^|)) and min 2 (maxt (M)), as 
functions of the propagation constant, \i. 



directions, cf. Ref. [4Tj : 

(^m,n H~ Wm,n) 



exp 



n ' &m,n 



XT] sech (rjt) , 
with r] taken as per Eq. (J3|, and 



( - -\ 
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a = In (2 (6 + m)) ■ 



V3 



-an, 



(4) 



(5) 



(6) 

a being the width of the hexagonal cell, see Fig. [l] 
The model's scale is fixed by setting a = 1. Factor 
(%m,n + Wm,n) m Eq. Q obviously corresponds to vor- 
ticity 5=1, and the exponential factor with a taken 
as per Eq. ([6| is determined as in the 2D spatial soli- 
ton with propagation constant — ja. The choice of 77 as 
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FIG. 3. (Color online) The phase and energy (f^™ \u m ,n\ 2 dt) 
patterns corresponding to input ansatz Q which generates an 
(unstable) hexagonal vortex with \i — 350. 



per Eq. (|3| implies that, simultaneously, the wave field 
in the factorized ansatz is localized in the longitudinal 
direction, in each core, as in the temporal soliton corre- 
sponding to the same propagation constant, cf. the struc- 
ture of the fundamental solitons considered above. The 
phase and energy structure of ansatz Q is illustrated, in 
a schematic form, by Fig. [3] 

Direct simulations of Eq. ([I]) with this input have been 
run in a broad range of values of the propagation con- 
stant, 7 < /i < 500. Nonetheless, stable spatiotempo- 
ral vortices with the hexagonal structure and topological 
charge 5 = 1 have never emerged. In fact, the evolution 
of the input organized as the ansatz of this type never 
leads to formation of any stable pattern. In the interval 
of 7 < \i < 13, the system makes an attempt to gener- 
ate a robust pattern of a triangular shape, as shown in 
Fig. [4j at three sites belonging to the original hexagon, 
the field quickly decays, while at three others it survives, 
for a while. However, the largest amplitude of the tem- 
porarily emerging triangular set is c± 8 (it is attained at 
fi = 13), while triangular vortices may be stable for am- 
plitudes above a threshold value of the amplitude which 
is ~ 18 (see below), therefore the triangles developing 
from the unstable hexagons are also subject to an insta- 
bility, eventually splitting into uncorrelated single-core 
excitations which separate in the longitudinal direction, 
see Fig. |4j 

Further, in the interval of 13 < fi < 70, the instability 
splits the original hexagon into a set of separating single- 
core excitations, the number of which varies randomly 
between 2 and 6 (not shown here in detail). In an adja- 
cent interval, 100 < fi < 200, the instability-development 
scenario is similar but faster, so that the formation of a 
transient triangular structure cannot be identified. 

At largest values of the propagation constant, 200 < 
ji < 500, the instability-development scenario is differ- 




FIG. 4. (Color online) The evolution of the six main compo- 
nents of the hexagon with initial propagation constant fi — 9. 
In this case, the simulations demonstrate the decay of the 
hexagon into a transient triangle, which is followed by a lon- 
gitudinal instability (splitting). 



ent. The six sites forming the hexagon keep their posi- 
tions and amplitudes for a while, but loose the mutual 
phase coherence. Then, instabilities of amplitudes and 
positions set in, but they manifest themselves on a much 
longer scale of the propagation distance, with z ranging 
from 10 to a few hundreds, instead of z ~ 1 at small /i, 
cf. Fig. |1J The separation between excitations in indi- 
vidual cores grows very slowly too, in comparison with 
the quick split of the transient triangle observed in Fig. 

H 



IV. THE GENERATION OF STABLE 
TRIANGULAR VORTICES 

The next step is an attempt to generate a triangu- 
lar vortical structure, which is suggested by the emer- 
gence of a transient one in the course of the evolution of 
the unstable hexagon (Fig. [4|. For this purpose, we 
used the same input as defined by Eqs. Q-^, but 
with three main peaks suppressed, which was done by 
replacing the fields at the corresponding sites by those 
from adjacent sites in the outer layer: w_i,o — > ^-2,0? 
u\^i — >• ^2,2? ^0,-1 ^0,-2? as shown in Fig. |5j The so 
constructed triangular ansatz keeps the vorticity of the 
original hexagon, 5=1. 

The evolution of this input was simulated in a broad 
range of values of the propagation constant, 4 < \i < 450. 
At small values, \i < 6, the three main peaks forming 
the triangle merge into a single-core fundamental soli- 
ton, which may be localized at the central site, or at 
any one belonging to the original triangle. Thus, stable 
vortices do not emerge in this case. For 7 < fi < 181, 
the triangle is destroyed by an instability which splits it 
into separating uncorrelated single-core excitations (the 
instability develops faster with the increase of /i). 

Finally, the same input generates stable triangular vor- 
tices at ji > 182, an example of which is shown in Fig. 
[6] The stability was verified by direct simulations for 
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FIG. 5. (Color online) The reduction of the hexagonal in- 
put to the triangular- vortex one: three of the six main peaks 
are replaced by fields taken from the surrounding layer, as 
shown by arrows (the replacement makes the amplitudes at 
the corresponding sites much smaller, but does not alter their 
phases) . 

long propagation distances, e.g., z = 986 for fi = 182. 
The temporal pulses in the cores which represent ver- 
tices of the triangle remain well phase-locked, keeping the 
phase circulation of 27r, which corresponds to vorticity 
5=1. The excitations at secondary sites (between the 
vertices) are phase- locked to the primary ones, but fea- 
turing some oscillations. The oscillations enhance with 
/i, but the overall vortical phase pattern always persists. 
On the other hand, the amplitudes of excitations at the 
secondary sites feature fast irregular oscillations, which 
also become stronger at larger \i (variations of these am- 
plitudes by a factor ~ 2 are observed already at the sta- 
bility threshold, /i = 182). These amplitude oscillations 
are coupled to small variations of amplitudes at the pri- 
mary sites, as shown in Fig. [7] It has been checked that 
the oscillations were not induced by reflection of small- 
amplitude radiation waves from edges of the integration 
domain (absorbers installed at the edges do not suppress 
the oscillations). 

Figure [8] displays the total energy, E = 
n f-tt \ u m,n(t)\ 2 dt, and amplitude of the tri- 
angular vortices, both stable and unstable ones, as 
functions of the effective propagation constant, \i + 
where the contribution S/i from oscillations of the 
fields is computed as follows. For each vertex of the 
triangle, (m,n), peak time £ m>n is defined, such that 
\um,n(tm,n)\ = max t |ix m>n (£) | , and the corresponding 
phase, 4>m,n(tm,n) ? is identified. Next, we compute 

where the average is taken over the three vertices of the 
triangle (or six ones for stable hexagonal vortices with 
5 = 2, see the next section), and, for the stable trian- 
gular modes, over z. For unstable triangles, the latter 
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FIG. 6. (Color online) (a) The transverse energy 
(f-tt \ u m,n\ 2 dt) and phase profile of the stable triangular vor- 
tex generated by the input with \i — 200. (b) The longitudinal 
(temporal) profile of excitations in the cores representing the 
vertices of the triangle. 

average was taken over a short interval Az, within which 
the pattern was not disturbed by the instability. As con- 
cerns the amplitude shown in Fig. [5J it was defined as 
\um,n (tm,n)\i averaged over z and over the three vertices, 
to smooth effects of small persistent oscillations of the lo- 
cal amplitudes. A small gap between the unstable and 
stable portions of the amplitude plot in Fig. ^ is due 
to the difference in wavenumber shift ([7]), as computed 
for the stable and unstable solutions. 



V. HEXAGONAL VORTICES WITH THE 
DOUBLE TOPOLOGICAL CHARGE 

As said above, the hexagonal input based on Eqs. Q- 
(|6| could not produce any stable pattern. However, the 
same initial ansatz, but with inverse signs of three of its 
main peaks — say u_i 5 o, ^1,1, uo-i — can give rise to sta- 
ble hexagonal spatiotemporal patterns carrying vorticity 
S = —2 (if the topological charge of the original ansatz is 
defined as S = +1), see an example of the stable mode in 
Fig. [9j Note that no changes were made to the hexagonal 
input at sites in outer layers. 

Simulations with this input were also run in a broad 
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FIG. 7. (Color online) The evolution of the longitudinal (tem- 
poral) profiles of one primary and one secondary components 
of a stable triangular spatiotemporal vortex, generated by the 
triangular input with \i — 450. 
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FIG. 8. (Color online) The amplitude and energy of the stable 
triangular vortex vs. the effective propagation constant, \i + 
see Eq. |7|. Blue (dashed) and red (solid) segments 
designate unstable and stable solution families, respectively. 



range of values of the propagation constant, 5 < \i < 450. 
At \i < 223, the seeded pattern is subject to various in- 
stabilities: spread out due to the lattice diffraction and 
temporal dispersion, or splitting into separating excita- 
tions, or, sometimes, merger into one or two single-core 
fundamental solitons. In particular, in the interval of 
8 < /i < 10, the pattern forms a transient triangular 
structure, which eventually splits, and in a broad inter- 
val of 20 < \i < 222 the initial hexagon fissions into two 
triangles, which also turn out to be unstable — essentially, 
because the amplitudes of the triangular patterns fall be- 
low the stability threshold. 

Stable double (S = 2) hexagonal vortices emerge at 
H > 223, whose shape is illustrated by Fig. [9] The 
stability of these vortices was verified in very long sim- 
ulations. With the further increase of /i, an instability 
island was revealed around \i = 300. In that case, the 
six temporal pulses remain locked to their positions, but 
the phase structure is lost at z > 500. Nevertheless, the 
hexagonal vortices recover their integrity at still larger \i. 
It is possible that other narrow intervals of the instability 
may be found inside the stability region. 



Figure 10 presents the energy and amplitude of the 
double vortices as functions of the effective wavenumber, 
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FIG. 9. (Color online) The same as in Fig. [6] but for a stable 
hexagonal vortex with topological charge S — 2, generated by 
the modified input with \i — 223. 
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FIG. 10. (Color online) The same as in Fig. [8] but for the 
hexagon-shaped vortices with the double topological charge. 
Only stable modes are presented in these plots. 



similar to Fig. [8] for the triangular vortices with 5 = 1. 
However, only the stable family of the hexagonal vortices 
is shown here, as we were not able to measure characteris- 
tics of unstable ones at fi < 223. In fact, the simulations 
produce no evidence that such unstable modes exist. 



VI. COMPACT TRIANGULAR VORTICES 

Still another type of stable spatiotemporal patterns 
can be produced by the input taken as per Eqs. Q-(|6|), 
but centered at an edge of the original hexagon, i.e., with 
m and n replaced by m — 2/3 and n — 1/3, respectively. In 
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FIG. 11. The same as in Fig. [5] but for a stable compact 
triangular vortex with topological charge S = 1, generated 
by the shifted input with ji = 15. 



this case, the results were collected for 5 < \i < 220, and, 
as shown in Fig. [TT] the stable structure takes the form 
of a vortex with 5 = 1, shaped as a densely packed trian- 
gle, without an empty site in the center, cf. Fig. [6] This 
structure is stable for /i > 11. It is relevant to stress that 
this stabilization threshold is more than an order of mag- 
nitude lower than its counterparts for the triangular and 
hexagonal spatiotemporal vortices reported in the previ- 
ous sections (recall those thresholds were /itriangie = 182 
and /ihexagon = 223, respectively). 

In fact, direct simulations initiated by the above- 
mentioned shifted input ansatz generate the compact tri- 
angle which seems "noisy". The noise can be removed 
by means of the "temporal filtering", setting the field 
equal to zero outside of the main pulse in each core, and 
running the additional propagation over Az = 10. Fur- 
thermore, for 31 < \i < 60, direct simulations starting 
from the shifted input ansatz lead to a phase instability. 
For instance, at \i = 32 and 40, the phases of the three 
vertices would take values 0, tt/2 and 7r, instead of those 
displayed in Fig. [6j Actually, this instability is caused by 
the fact that the input is far from the shape of the sta- 
ble mode, giving rise to several temporal peaks in each 
core. If the initial data are "cleaned up" by nullifying 
the field outside of the main temporal pulse, the simula- 
tions converge to stable compact triangular vortices. The 
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FIG. 12. Trajectories of the motion of colliding triangular 
vortices rotated by angle tv/3 relative to each other, for \i — 
250 and velocities ±ko = ±11. 



amplitude of stable compact triangular vortices evolves 
slowly and almost linearly versus the effective propaga- 
tion constant fi + Sfi: the computed values of the latter 
range from ~ 22 to 78, then the amplitude goes from 40.3 
to 41.2. 



VII. COLLISIONS BETWEEN MOVING 
VORTEX SOLITONS 

The availability of stable solutions for the vortex spa- 
tiotemporal solitons, and the obvious Galilean invariance 
of Eq. ([I]) suggest to study collisions between moving 
vortices. In particular, it is interesting to simulate colli- 
sions between stable triangular modes shown in Fig. [6j 
rotated by angle 7r/3 relative to each other, to test a pos- 
sibility of their fusion into a full hexagonal vortex of the 
type displayed in Fig. [9] This was done taking a pair 
of the triangles separated by a relatively large temporal 
interval, At = 16, for values of ji > 182, at which the 
triangular vortices are stable by themselves, as shown 
above. They were set in motion, multiplying them by 
exp(±ifco£), which, obviously, lends the solitons veloci- 
ties ±ko (in terms of the optical waveguides, these are 
shifts of the inverse velocities). 

In fact, the fusion of colliding triangles into a hexagon 
was never observed. Instead, slowly moving triangles 
demonstrate a long-range repulsion and stop at finite dis- 
tance (but do not bounce back), as shown in Fig. 12 At 
intermediate velocities, the colliding triangular vortices 
do bounce back, and eventually they get destroyed by the 
longitudinal instability (splitting into uncorrelated tem- 
poral pulses in different cores), as shown in Fig. [l3j At 
high velocities, the solitons, quite naturally, pass through 
each other, loosing some kinetic energy. There is a sharp 
threshold between the rebound regime and the passage. 
Just above this threshold, the passing vortices get de- 
stroyed by the longitudinal instability shortly after the 
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FIG. 13. The collision of two triangular vortices at \i — 250 
and ko = 12. Trajectories of individual pulses forming the 
vortices are displayed. 
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remain stable indefinitely long after such quasi-elastic 
collisions. Indeed, since the numerical box has a finite 
length, and periodic boundary conditions in t are used, 
the moving vortices undergo repeated collisions, loos- 
ing some velocity each time. Eventually, they would be 
destabilized by a collision occurring at a lower speed. 

VIII. CONCLUSION 

We have introduced a system of parallel waveguides 
with the linear coupling between nearest neighbors, based 
on the hexagonal lattice in the transverse plane. Each 
guiding core features the cubic self-attractive nonlinear- 
ity. The system can find straightforward realizations 
in nonlinear optics, and in BEC trapped in the corre- 
sponding optical lattice. Systematic simulations, start- 
ing with a natural input ansatz for vortical hexagons, 
reveal three distinct species of stable semi-discrete spa- 
tiotemporal complexes, which are discrete in the trans- 
verse plane and continuous in the longitudinal direction. 
These are triangular modes with vorticity S = 1 and 
hexagonal ones with 5 = 2, both built with an empty 
core at the center, and compact triangles carrying 5 = 1, 
without the central empty core. Collisions between stable 
triangular vortices were also studied by means of simula- 
tions, demonstrating the stoppage of the slowly moving 
vortex solitons, destabilizing rebounds, and quasi-elastic 
passage, depending on the collision velocity. 

More complex structure of the arrayed waveguides can 
be considered as a generalization of this work (in partic- 
ular, quasi-periodic lattices). It may also be interesting 
to study vortex complexes in two-component models of 
the same type. 



FIG. 14. (Color online) Regions of different outcomes of col- 
lisions between two mutually symmetric triangular vortices, 
rotated by angle ty/3. 



collision. Domains corresponding to different outcomes 
of the collisions in the (/i, ko) plane are shown in Fig. 14 
While the collisions seem elastic with the increase of 
the velocity, it was not possible to conclude if the vortices 
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